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1. Introduction 



Consider an inextensible string of unit tension placed along the x axis with its left 
endpoint at 0, and let M(x) be the total mass of the string segment from to x. The 
small vertical oscillations y = y(x) of the string obey the equation 

y" + zM'y = (1.1) 

where z is the square of the oscillation frequency. The natural frequencies of the string 
depend on how it is tied at the ends; suppose 

y'(O-) = y(£) + (L - £) y'(£+) = (1.2) 

where i is the supremum of the points of increase of M, and i < L < oo is some "tying 
constant" . 

Now, let z G R_, and denote by (f(-,z) and ip(-,z) the two particular solutions of 
Equation (jl.ip satisfying 



(p(0 t z) = V>'(0- z) = 1 and (f/(0-, z) = ip(0, z) = 
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The characteristic function of the string is then defined by 

W(z) = lim ^4 (1-3) 
or, equivalent ly, 

WW = (1-4) 
w'(0—, z) 

where w(-,z) is, up to an unimportant factor, the unique non-negative solution of the 
string equation that is decreasing and satisfies w(L, z) = 0. The definition of W may be 
extended to C\M + by analytic continuation; W is the analog of the Weyl-Titchmarsh 
function in the Sturm-Liouville theory. 

The present paper is devoted to the inverse spectral problem for this vibrating string 
equation, namely the problem of recovering M from W. We call this Krein's inverse 
problem. We describe and illustrate a novel algorithm that is applicable to a class of 
characteristic functions associated with a certain continued fraction. The remainder 
of this introduction provides a summary of the approach and discusses its relationship 
with other works. 

1.1. Krein's inverse problem 

In a series of papers published in the Soviet Union during the 1950's, M. G. Krein made 
a detailed study of the existence of spectral expansions associated with the vibrating 
string equation. For some useful accounts of this work in the English language, see 
[T0| IT5| fl8] : in particular, [15] contains a rough outline of the historical development 
of Krein's ideas and their overlap with the work of Feller aimed at a unified analytical 
treatment of some discrete and continuous stochastic processes. 

Denote by M the set of functions M : [0, oo] — > [0, oo] that are non-decreasing, 
right-continuous and infinite at infinity. We use the convention M(0— ) = so that 
dM(x) is a well-defined Stieltjes measure. The tying constant L is absorbed in the 
definition of M by setting 

L := sup {x E R+ : M(x) < oo} . 

For such M, Krein gave a precise meaning to Equation (11.11) and its particular solutions 
ip(-, z), ip(-,z) and w(-,z). With a slight abuse of terminology, we shall identify M with 
the string. Krein showed that the characteristic function defined by equation (11.31) is 
necessarily of the form 

f d °W s f°° d(r W 
w{z) = c+ -v^j = c + — + y-4 ( ) 

./[0,oo) A ~ Z ~ Z JO A — Z 

where c > 0, a is a right-continuous non-decreasing function with support in [0, oo) such 
that 

TTT < - M 

[0,oo) 1 + A 
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and s = cr(0). The paramater c is in fact the infimum of the points of increase of M, 
and the function o in this representation is called the principal spectral function of the 
string. We shall denote by W the set of all functions expressible in this form. Hence 
Krein showed that W G W for every M G M. He also conjectured the converse: every 
element of W is the characteristic function of a unique string in M. To the best of our 
knowledge, the first proof of this conjecture was given in [10], Chapter 6; it uses the 
theory of Hilbert spaces of entire functions, and a deep uniqueness theorem due to de 
Branges. 



1.2. The Sturm-Liouville problem in impedance form 

Although there is no known systematic procedure for solving Krein's inverse problem 
in its full generality, effective algorithms do exist if the mass distribution belongs to 
certain special classes. For instance, if M is smooth, then the change of variables 



— [q— ) + zgu = (1.7) 



t = T{x):= / y/M'(s)ds, u{t) = y{X{t)) 
Jo 

where X is the inverse of T, transforms the string equation into 

d ( du\ 

dt\ ~dt) 
where 

e = VM> o X (1.8) 

is the so-called impedance. The inverse problem for this equation may be tackled by 
"integral operator" techniques going back to the works of Gelfand, Krein, Levitan and 
Marchenko; see [19] for an informal account of these techniques and their historical 
background, and [1] for more recent developments. General methods inspired by the 
Gelfand-Levitan approach aimed at solving the inverse problem for the string equation 
have been considered by Dym & Kravitsky [8j|9], Winkler [23] and Boumenir [6]. 

1.3. The Stieltjes moment problem and the class M + 

Another class for which the inverse problem is tractable arises in connection with 
Stieltjes' moment problem [20J: given a sequence of numbers Cq, c±, . . ., find a measure 
da on R + such that 

POO 

c n = X n da(X) for n = 0,1,.... (1.9) 
Jo 

Briefly, Stieltjes' approach was to construct from the Cj a certain continued fraction, 
namely 

l — • (i.io) 

s z H j 

si + 

-s 2 z H 
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Stieltjes gave an algorithm that permits the calculation of the s n from the c n . The 
coefficients s n thus obtained form a finite or infinite sequence of strictly positive numbers, 
depending on whether the measure a has finitely or infinitely many points of growth. 
Stieltjes distinguished the determinate case, in which the moment problem has a unique 
solution, and its complement, the indeterminate case. The determinate case arises if 
and only if the continued fraction either terminates (i.e. the sequence of the s n is finite) 
or else is infinite but convergent. The continued fraction is then a well-defined function 
of z which may be expressed as the Stieltjes transform of the sought measure a . Stieltjes 
also found a beautifully simple criterion for the convergence of the continued fraction 
when the sequence of the s n is infinite; he showed that convergence occurs if and only if 

oo 

^s n = oo. (1.11) 

n=0 

This work on the moment problem provided the inspiration for Krein's own theory 
of strings; see Supplement II in [H]. With remarkable insight, Krein was able to relate 
Stieltjes' continued fraction to the characteristic function of a discrete string, i.e. a 
string made up of a sequence of point masses rrij at the positions Xj. To see how, set 

m/ \ W(X,Z) 

W{x,z) := -. 

w'{x, Z) 

The string equation then yields a Riccati equation for W(-,z): 

W'(x, z) = -l- zM'(x)W 2 (x, z) . 

For a string consisting of point masses, M' is a linear combination of Dirac deltas; at a 
point Xj with concentrated mass mj, this Riccati equation should be interpreted as 

W(xj—,z) = W(xj+,z) + zmjW(xj-,z)W(xj+,z) . 

By using Equation (11.41) . we find 

W(z) = W(0-, z) = x + W(x -, z)=x + = 



-m z + 



W(x +,z) 



and, by iterating, 



W(z)=x + j . (1.12) 

-m z H 



Xl - Xq H 1 

-m n z + 



W(x n +,z) 

If the principal spectral function has finitely many points of increase, then the string 
consists of finitely many, say n, point masses with finite mass (and a further point of 
infinite mass at L). Then W(x n +, z) = L — £, the expansion terminates, and we have 
the finite case of Stieltjes' continued fraction (jl.lOp . If the principal spectral function 
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has infinitely many points of increase, then there are infinitely many point masses; the 
continued fraction does not terminate but its coefficients obviously still determine M(x) 
for every x such that 

x < £ := lim x n . 

In particular, if 

oo 

£o + TTlj = OO 

then the continued fraction converges and the string is completely determined by its 
coefficients since either £$ = oo or else L = £q. 

Let us now elaborate the significance of the foregoing remarks for Krein's inverse 
problem. Define 

jVeW:^ \ n do-(\) <oo for n = 0, 1, 2, ...J (1.13) 

and denote by M + the corresponding set of strings. Every W G W + may be expanded 
in a continued fraction of the form f ll . 12j) . and the coefficients in the expansion may 
be computed from the moments of the principal spectral function o. It will serve our 
purpose to abuse Stieltjes' terminology somewhat by saying that W itself belongs to 
the determinate subclass of W + if the continued fraction is either finite or convergent, 
and to the indeterminate subclass otherwise. It follows, then, that Stieltjes' continued 
fraction provides an effective solution method of the inverse problem for strings in M + 
whose characteristic functions are in the determinate subclass. 

Remark 1. Not every member of M + is a discrete string; all that can be said is that 
such strings must begin with a sequence of point masses; see [7^ /. §5.9. If W G W + 
is in the determinate subclass, then M is a discrete string. The converse is of course 
not true, and for a discrete string whose characteristic function is in the indeterminate 
subclass, one cannot recover the tying constant from the continued fraction coefficients 
alone; see Uhl. §13. 



l.Jf.. Main results and outline of the paper 

The connection between discrete strings and continued fractions can be exploited to 
provide a reconstruction algorithm of greater applicability. Beals et al use this idea in 
the context of integrable systems (21 E] ; Borcea et al use it as the basis of their numerical 
treatment of the inverse Sturm-Liouville problem in impedance form [4j. The practical 
issue that arises is how to find approximations of the spectral data that satisfy Stieltjes' 
moment condition and have explicitly computable moments. Both Beals et al and 
Borcea et al resolve this issue by restricting their attention to the case of a discrete 
spectrum, where the characteristic function has a Mittag-Leffier expansion which may 
be truncated to furnish the required approximations. 

We present an alternative truncation strategy applicable to a class that includes 
some characteristic functions with an absolutely continuous spectrum. This class is 
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a proper subset of a class denoted M~, based on a continued fraction expansion that 
arises naturally in connection with Krein's inverse problem. The expansion rests on the 
following two useful "rules" relating the characteristic function to its string. The first 
of these rules concerns the characteristic function of the right-continuous inverse M* of 
the string M: 



Proposition 1.1. 



W*(z) 



-zW(z) 



The string M* will be called the dual of M. The second rule concerns the string 
obtained from M by removing from the spectrum of a: 

Proposition 1.2. Let M be the string with characteristic function 

da(A) 



Then 



W(z) 



W{z) :=c + 



c + 



[0,oo) 



X-z 



d " (A) - W(z) + S - 



X-z 



is the characteristic function of the string M defined by 



M(x) :-- 



M{t) 

M(t) 
M(oo-) 



where x 



M(t) 
M(oo- 



n 2 



dr. 



Two other useful formulae connecting W and M are: 

M(oo-) = - (1.14) 
s 

and, in the case t + M(£—) < oo, 

L = c+ [ d ^>. (1 . 15) 

J[0,oa) X 

For proofs of these results, see [TO] , §6.8 and §6.9. Using these rules and setting 
Wo(z) := W(z), we may write 

W (z) = ^ + W (z) = ^ + 



-z uw -z -zW^z) 

where W\ is the characteristic function of the dual string, say Mi, of M . By the same 
argument, 

. Sn 1 
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where W2 is the characteristic function of the dual string, say M 2 , of Mj. Iterating, we 
obtain 

W (z) = ^ + j . (1.16) 

Z si + • • • + 1 



-zW 2 „-i(z) 

This continued fraction, albeit in a very particular context and in a somewhat disguised 
form, was introduced and studied in our recent work on the excursions of diffusion 
processes [TJ. It is also implicit in some of the calculations carried out by Donati- 
Martin and Yor in their generalisation of Levy's formula for the area enclosed by a 
planar Brownian motion [11] . 

In £j2j we study the relationship between this continued fraction and that of Stieltjes. 
This relationship may be described in terms of the map 

W ^ W~(z) := —W* (-) . (1.17) 



Then 



-z 



W(z) = — + = ^ W~(z) = 

— z 1 1 
si + s z + 

—z 



Sl + 



-s 2 z 



Hence the class 

W" := {W~ : W G W + } (1.18) 

is the set of characteristic functions that plays, for the continued fraction (I1.16p . the part 
played by the class W + for the continued fraction of Stieltjes. We show by elementary 
means that the corresponding set of strings, denoted M~, consists of those strings whose 
dual principal spectral functions have negative moments. Then, by a straightforward 
application of Kats' results on such spectral functions [TB], we derive a simple criterion, 
expressed in terms of the mass distribution M(x), for M to be in the class M~. The 
upshot is that this class is large enough to be of interest; in particular, it contains some 
strings with an absolutely continuous principal spectral function. 

Just as in the case of the Stieltjes continued fraction, W~ may be partitioned into 
two subclasses: determinate and indeterminate. The determinate subclass consists of all 
those characteristic functions for which the continued fraction expansion (I1.16P is either 
finite or else convergent — the latter case occuring if and only if Equation ( 11. lip holds. 
§|3] discusses the finite case. We show there that the string M corresponding to such 
a characteristic function is discrete with finitely many point masses, and we devise an 
algorithm to compute it explicitly, given the continued fraction coefficients. The infinite 
(convergent) case may then be tackled by considering the truncations of the continued 
fraction after n terms; since Krein's correspondence between strings and characteristic 
functions is a homeomorphism, the discrete string so obtained yields in the limit n — > 00 
the solution of the inverse problem. 
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In §U we give some examples that arise in the study of the excursions of one- 
dimensional diffusion processes. Knight |T7] and Kotani & Watanabe [18] pointed out 
simultaneously the relevance of Krein's theory to this topic: every string M defines 
a generalised diffusion process, and the excursions of the process from its starting 
point may be described by means of the principal sprectral function of the dual string. 
The probabilistic version of Krein's inverse problem is to find the diffusion, given the 
distribution of the excursion lengths. Donati-Martin and Yor solved this problem 
explicitly in a number of interesting cases (T2J 13] , and we use their findings to illustrate 
the effectiveness of our algorithm. 

The paper ends with a few concluding remarks in £j5j 

2. The class M~ 

Lemma 2.1. For every W G W, W~ G W. Furthermore, the map 



W 



from W to itself is an involution. 



Proof. Let W G W. Since W* G W, we can write 




where c* > 0, s* > and 




A simple calculation shows that 




where 



da~(\) 



da*(l/\) 
A 



(2.1) 



It is then readily seen that 




and so W~ G W. The fact that 



(w-) 



w 



then follows easily. 



□ 



Theorem 1. 



W 
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Proof. The fact that W h- > H 7- is an involution implies 

W" := {W : W G W + } = [W G W : W" G W + } 
Note that, by definition of W + , W~ G W+ if and only if 

A" da' (A) < oo for n — 0,1,... 



The theorem is then a consequence of Equation ( 12. ip . □ 
Corollary 2.1. 

MeM^ ^ M* g Mr . 



Proof. The dual of a discrete string is a discrete string. Hence 

M G M + <^ M* G M + . 

Furthermore, the maps 

M v+ M~ and M H- M* 
commute. □ 

Next, we seek a simple criterion — in terms of the mass distribution — for M to 
belong to the class M~. We shall need the following key result, due to Kats [T6] : 

Kats' Theorem. For the principal spectral function a of a string M G M to satisfy 

/ A~ n dcr(X) < oo for every n = 1, 2, . . . 

</[0,oo) 

it is necessary and sufficient that L < oo and that 

V e > 0, M(x) [L - x} 1+£ = o(l) asx^ L-. 

Theorem 2. M G M~ i/ and on/?/ i/ either 

(i) M(oo-) < oo and We > 0, linr,^ x [M(oo-) - M(x)] 1+£ = 
or 

(ii) L < oo and V £ > 0, lim x ^ L _ M(x) [L - x] 1+£ = 0. 

Proof. We note first that, for M G M~, the implication 

s = s* > 

holds because W _ has an expansion of the form (j!.12p . Furthermore M(oo— ) = K 
Hence we have the trivial identity 

M G MT [M G MT and M(oo-) < oo] or [M G Mr and M(oo-) = oo] 



Krein's inverse problem and continued fractions 10 

[M G M~ and M(oo-) < oo] or [M* G M~ and M*(oo-) < oo] . 
Let us show first that 

[M* G M~ and M*(oo-) < oo] Condition (ii) holds. (2.2) 

Suppose that M* G MT and M*(oo-) < oo. Then 

£ = L = M*(oo-) < oo 

and so 

s = 0. 

It follows from Corollary 12.11 that 



/ \- n da{\) < oo forn = l, 2, .... 

ifO.oo) 



' [0, oo) 

Kats' Theorem then implies that 

Ve>0, M(x)[L-x] l+£ = o(l) asx-^L-. 

This is Condition (ii). Conversely, if Condition (ii) holds, we arrive by a similar argument 
at the conclusion that M* G M~ and M*(oo— ) < oo. This proves (12.21) . To complete 
the proof of the theorem, it suffices to replace M by M* in f)2.2p . □ 

Remark 2. The problem of characterising the determinate subclass in terms of the 
mass distribution appears to be much more difficult. 



Example 1. Let «(■, z) be the particular solution of Equation (1.1) in the interval [0, 1] 
such that 

u'(0, z) = —1 and u(l, z) = . 

The problem considered by Borcea et al is to recover g, normalised by g(0) = 1, from 
u(0, z) W. Now, suppose that 

1 dr f 1 

< oo and / g(r) dr < oo . 



/o e( T ) jo 

Then u(0, z) is the characteristic function of a string M such that 

i = L < oo and M(L-) < oo . 

It is an easy consequence of Theorem IE that this string belongs to M~ . Knowing M, g 
can in principle be obtained via Equation U.8\) . 

Given the continued fraction coefficients s n , one can ascertain whether W is in the 
determinate class by using Stieltjes criterion M.ll\) . To give a more concrete example, 
consider the case 

0=1. 
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Then 

W(z) = ; = 1 . 



1 + 



3 1 
+ 



-z 5 H 

The series 

oo oo 

5> = £(2 J + 1) 

n=0 j=0 

diverges, and so W is in the determinate class. 
3. The inversion algorithm 

We proceed to compute the string corresponding to the finite case of the continued 
fraction expansion on the right-hand side of Equation (11.161) . In what follows, it will 
always be assumed that 

s > and sj > if j > . 
We shall use the compact notation 



W[Sj, . . . , S j+2k ] ■= — 



: 1 



J + i S,- +2 1 



1 

Sj+2k-l + 

— Z 



and 

W[Sj, . . . , s j+2k +i] ■= — + 



z 

Sj+1 



—Z 



1 

s j+2fc-i + ; f- 

s i+2fc 1 



and an analogous notation for the corresponding strings. These characteristic functions 
are in the determinate subclass of W~ because their images under the map W i— > W~ 
are in the determinate subclass of W + . 

Theorem 3. For every k > 0, 

M[s k ] = - 

and, for every < j < k, 
where 

x= I {1 + Sj M*[s j+1 , s k ]{r)} 2 dr. (3.2) 
Jo 
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Proof. Consider the case where k — j is odd. Then 

W[ Sj , ...,s k } = ^ + — . 

+ 1 r 

Sj +2 i 



— + ■■■ + 

— z 



-z s k 



An easy calculation shows that 

(w) [s j} ■ ■ ■ , s k ] = W[s j+U ...,s k ]. 
By taking the dual on both sides, we obtain 

W[sj, . .., s k ) = W*[s j+1 , s k ]. 
It then follows from Proposition 11.21 that 

M[ Sj , s k ](t) 



M*[s j+1 , . . . , s k ](x) 



1 - SjM[sj, . . . , s k )(t) 
where 



x 



[ {l- Sj M[ Sj , . . . , s fc ](r)} 2 dr . 
Jo 



When we "turn this around" and express M[sj, . . . , s k ] in terms of M*[sj+\, . . . , s k ] 
and t in terms of x, we obtain the desired result. The case where k — j is even is 
analogous. □ 

Corollary 3.1. For every k > and every < j < k, M[sj, . . . , s k ] is a discrete string. 

Proof. Let k > and proceed by induction on j, starting with j = k — 1 and going 
backwards. The claim is obviously true for j = k — 1. Suppose that M[sj + x, . . . , s k ] 
is a discrete string. By Equation (13. 2p . x is then an increasing continuous piecewise 
linear function of x, and the fact that M[sj, . . . , s k } is piecewise constant is an obvious 
consequence of Equation (13. ip . □ 

This result permits the efficient recovery of the discrete string M[sq, . . . , s n ] via 
the sequence 

M[s n ] -> M[ ] ->■ ► Af [So, • • ■ , s n ] . (3.3) 

Each term in that sequence is a right-continuous piecewise constant function with finitely 
many jumps, so we need only keep track of the location, say Xj, of these jumps, and of 
the value, say yj, of the function there. We always include x = and y , whether or 
not y > 0. We have 



So we set 



M[s n ](x) = — for every x > . 
xf = and yf = — . (3.4) 

S m. 
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M[s n _i,s„](x) = 







if < x < - 1 - 
if x > -i- 



Sn-l 



So we set 



x, 



(i) 



0, ac^ = — and y£> = 0, y] 1 



.(!) 



.(!) 



Then, for > 1, we have the recurrence relations 



x, 



(2k) 



X 



(2k) _ (2k) 



x)_{ + 



1 + S n -2k Xj 



(2k-l) 



(2k-l) (2k-l) 

y) -v)-i 



(2k) 



X 



(2k-l) 
j+1 



1 ~r 6 n -2fe x j+i 



0<j<k. 



S n -2k 



X 



(2k+l) 




n (2k+l) (2k) 

0, x\ -Vo 



(2k+l) _ (2k+l) 

x j+1 - xj -r 



1 + Sn-2k-l x 



(2k) 



y 



(2k) _ (2k) 

y 7— 



?/ ( 2fe + 1 ) 



0. 



y 



(2k+l) 



X 



(2k) 



1 + S n _2fc-1 Xj 
1 



(2k) 



1 < j < k . 



l<J<k 
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(3.5) 

(3.6) 
(3.7) 

(3.8) 



(3.9) 
(3.10) 

l<j<k. (3.11) 
(3.12) 

(3.13) 



(2k+l) 

Vk+i = • 

The computation terminates when the superscript reaches the value n. 



(3.14) 



4. Diffusion processes 

In this section, we describe in broad terms the probabilistic version of Krein's inverse 
problem expounded by Knight in [17], and illustrate our inversion algorithm by means 
of an example taken from [12J. The reader unfamiliar with diffusion processes will find 
a useful summary of the theory in [3]. 

Consider a one- dimensional diffusion process X started at the origin, with values 
in some interval / C R + . Suppose that the origin is instantaneously reflecting and 
that X is in natural scale. Such a process is completely determined by its infinitesimal 
generator, i.e. by the differential operator 

1 d 2 



M'(x) dx 2 



(4.1) 
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acting on a suitable set of functions. Here, dM is the speed measure, and the domain of £f 
consists of twice differentiable functions /:/—»■ M + satisfying the condition /'(0+) = 0; 
an additional condition may be imposed on / at the right boundary, depending on the 
behaviour of the diffusion there. The process 

is called the local time (at the origin); it measures the time spent by the diffusion X in 
the vicinity of the origin, up to time t. Its right-continuous inverse 

r t : = inf {, : Inn — L_ ^ {Xu ) du > t 

is called the inverse local time (at the origin); it is a positive, non- decreasing process 
that jumps at the times when X begins an excursion away from the origin. The height 
of the jump is the length of the excursion, i.e. the time that elapses before X returns to 
its starting point. It may be shown that the inverse local time is in fact a Levy process. 
From the theory of such processes, one deduces the existence of a Levy exponent 0(A) 
defined implicitly by 

E (-An) = exp [-te(A)] for A > . (4.2) 

Furthermore, the Levy exponent is necessarily of the form 

poo 

0{X) = a\ + b+ (l - e~ Xy ) dv{y) (4.3) 
Jo 

for some numbers a, b > and some Levy measure du, i.e. a measure such that 



roc 

/ min{l, y} dv(y) < oo . 
Jo 



Knight shows that, if X is a so-called "gap diffusion", then the function M appearing 
in Equation (14. ip may be identified uniquely with a string, still denoted M. The Levy 
exponent is related to the characteristic function W of the string via 

Furthermore 

;>oo 

a = c *, b = s* and v\y) = / Ce~ C W(C) . (4.5) 

Jo 

Krein's inverse problem may thus be rephrased as: "Find the gap diffusion, given the 
Levy exponent of the inverse local time" . 

To give an example, let the Levy exponent be given by Equation (14.31) with 
a = b = and 

"'(y)--=C^-, «e(0,l), (3>0. (4.6) 

y l+a 
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Donati-Martin and Yor showed that the corresponding diffusion is, up to a 
homeomorphism, a Bessel process with drift [121 122]. The corresponding string may 
be expressed in terms of the modified Bessel functions. Set 

1 a/7 



W(z) :-- 



(4.7) 



e(- z ) {i-z/py-i 

where 

7 :=CT(l-a) P a . 

By using the well-known continued fraction expansion for the binomial (see for instance 
p. 343), it is easy to see that W is expressible in the form (jl.lfip with 



s 2 j - - 



a) 



- (2j + 1) and s 2j+ i = 27 



:i + « 



a 



where (•)„ denotes Pochhammer's symbol and we use the convention (•)„ = 1 for 
n < 0. Stieltjes' criterion (11. lip implies that this continued fraction converges. Hence 
W belongs to the determinate subclass of W~. In what follows, we find approximations 
of the string M by truncating the continued fraction after n terms, and computing the 
corresponding discrete string by the algorithm of the previous section. 

Example 2. Take 

' = 2 and C 



a 



1 

2' 



1 



'2tt 



Then 



M(x) 



2x 



(4.8) 



1 + 4x 

and the diffusion is a Brownian motion with drift —2 in natural scale JE/. The dots in 
FigureU\ are the points (xj,yj n ^) such that x™ < 5 corresponding to the discrete string 
M[sq, . . . , s n ] constructed by the algorithm of 53 The superimposed continuous curve 
is a plot of M. 



Example 3. By taking 



and letting a — » 0+ ; we obtain 

2 



C 



2Y{l-a)l3 a 



W(z) 



(4.9) 



\n(l-z/P) 

The string with this characteristic function was found explicitly by Donati-Martin 

and Yor JWj : let Y be the R + -valued diffusion process, reflected at the origin, with 
infinitesimal generator 



Y 



2^ 



+ 



1 

2^ 



K (V2j3y) 



K (VWy) 



d 

dy 



The scale function of Y is 



dt 
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0.5 n 




1 2 3 4 5 



X 

Figure 1. The points (af„ , y„ ), shown as dots, corresponding to the approximation 
M[sq, . . . , s n ] of the string M defined by Equation (|4.8p . Here, n = 511 and the solid 
curve is a plot of the true M. 



The diffusion 

X = S{Y) 

is in natural scale and its generator is given by Equation j4-l\) , where M is the string 
corresponding to Figured shows a plot of the string for x small, together with an 

approximation obtained by our algorithm. 

5. Concluding remarks 

We end with some comments on the numerical issues arising from the algorithm: by 
design, the approximation is exact in the limit x — > oo. More generally, our computations 
suggest that 

max \y$ - M (x^) \ = O (n~ 1/2 ) as n -> oo . 

There is also some numerical evidence that greater accuracy may be obtained by 
averaging at the jumps, namely 



max 

3 



yi ^ — M (x$) =0(n" 1 ) as n 



— > oo 



where 

(3-2) _ y n + v n 
y n 2 

As one would expect, in order to compute the discrete string from the truncated 
continued fraction, it is necessary to have good approximations of the coefficients. It 
will seldom be the case that exact formulae are available. Instead, one will need to 
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0.12- 
0.1- 
0.08- 

M 

0.06- 
0.04- 
0.02- 



u 0.2 0.4 0.6 0.8 1 

X 

Figure 2. The points {%n ,Un ), shown as dots, corresponding to the approximation 
M[sq, . . . , s n ] of the string M whose characteristic function is given by Equation (|4.9p . 
/3 = 2, n = 2047 and the solid curve is a plot of the true M. 

compute these coefficients by using some numerical algorithm; see for instance [7] and 
the comments therein. 

Finally, the inversion method we have described generates only piecewise constant 
approximations of M, and so the recovery of M' — in cases where the string is absolutely 
continuous — is not entirely straightforward. It would be of interest to adapt to our 
case the ingenious recovery technique that Borcea et al jl] devised in the context of the 
Sturm-Liouville problem in impedance form. 
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